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ABSTRACT 

Context. Massive stars form in dense and massive molecular cores. The exact formation mechanism is unclear, but it 
is possible that some massive stars are formed by processes similar to those that produce the low-mass stars, with 
accretion/ejection phenomena occurring at some point of the evolution of the protostar. This picture seems to be 
supported by the detection of a collimated stellar wind emanating from the massive protostar IRAS 16547—4247. A 
triple radio source is associated with the protostar: a compact core and two radio lobes. The emission of the southern 
lobe is clearly non-thermal. Such emission is interpreted as synchrotron radiation produced by relativistic electrons 
locally accelerated at the termination point of a thermal jet. Since the ambient medium is determined by the properties 
of the molecular cloud in which the whole system is embedded, we can expect high densities of particles and infrared 
photons. Because of the confirmed presence of relativistic electrons, inverse Compton and relativistic Bremsstrahlung 
interactions are unavoidable. 

Aims. We aim at making quantitative predictions of the spectral energy distribution of the non-thermal spots generated 
by massive young stellar objects, with emphasis on the particular case of IRAS 16547—4247. 

Methods. We study the high-energy emission generated by the relativistic electrons that produce the non-thermal radio 
source in IRAS 16547—4247. We also study the result of proton acceleration at the terminal shock of the thermal jet 
and make estimates of the secondary gamma-rays and electron-positron pairs produced by pion decay. 
Results. We present spectral energy distributions for the southern lobe of IRAS 16547—4247, for a variety of conditions. 
We show that high-energy emission might be detectable from this object in the gamma-ray domain. The source may 
also be detectable at X-rays through long exposures with current X-ray instruments. 

Conclusions. Gamma-ray telescopes like GLAST, and even ground-based Cherenkov arrays of new generation can be 
used to study non-thermal processes occurring during the formation of massive stars. 
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1. Introduction 

The mechanism of formation of massive stars remains one 
of the open questions in the held of star formation. It is 
known that these stars originate inside massive molecular 
clouds but the sequence of processes that take place dur- 
ing the formation of the star are mostly unknown. It has 
been suggested, for example, that the coalescence of vari- 
ous protostars in the same cloud can lead to the emergence 
of a massive star (e.g. Bonnell, Bate, & Zinnecker 1998). 
Massive stars appear in massive stellar associations where 
cloud fragmentation seems to be common. Alternatively, a 
massive star could form by the collapse of the core of a mas- 
sive cloud, with associated episodes of mass accretion and 
ejection, as observed in low-mass stars (e.g. Shu, Adams & 
Lizano 1987). In such a case, the effects of jets propagating 
through the medium that surrounds the protostar should 
be detectable. 
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Recently, Garay et al. (2003) have detected a triple 
radio continuum source associated with the protostar 
IRAS 16547—4247. The radio source presents a linear struc- 
ture consisting of a thermal core, and two radio lobes. The 
southern lobe is clearly non-thermal, indicating the pres- 
ence of relativistic electrons that produce the observed radi- 
ation by synchrotron mechanism. This non-thermal source 
has been interpreted by Garay et al. (2003) as the termi- 
nation point of one of the jets ejected by the protostar. 
There, a strong shock would accelerate the electrons up to 
relativistic energies by Fermi mechanism (e.g. Bell 1978). 
The observed spectral index of a ~ —0.6 (S u oc v a ) is in 
good agreement with what is expected from an uncooled 
population of relativistic electrons produced by diffusive 
shock acceleration at a strong non-relativistic shock (e.g. 
Protheroe 1999). 

The angular separation of the non-thermal source from 
the core corresponds to a linear distance of only 0.14 pc, 
and then the population of relativistic particles is inside the 
molecular cloud. Hence, these particles are in a rich envi- 
ronment, with a high density of ambient matter and pho- 
ton held from the infrared emission of the cloud. Inverse 
Compton (IC) and relativistic Bremsstrahlung losses are 
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then unavoidable for these particles. If protons are accel- 
erated at the termination shock along with the electrons, 
then inelastic pp collisions can take place, producing pions, 
which will decay yielding gamma-rays, relativistic electron- 
positron pairs and neutrinos. The radiation produced via 
all these mechanisms will be likely steady at scales of years 
due to the dynamical timescales of the processes occurring 
at the source. 

The main goal of the present paper is to estimate the 
high-energy yield of all these interactions, both leptonic 
and hadronic, in order to ponder whether gamma-ray as- 
tronomy can be used to probe the massive star formation 
and the outflows it could produce. Till now, thermal radio 
and X-ray emission has been associated with the formation 
of low-mass stars. Here we will show that massive proto- 
stars can produce a significant amount of radiation in the 
gamma-ray domain, because of the dense and rich medium 
in which they are formed. 

The structure of the paper is as follows. In the 
next section we provide the basic information about 
IRAS 16547—4247, the associated radio sources, and the 
ambient medium. Then, in Section [3] we discuss the parti- 
cle acceleration and the different losses for the relativistic 
particles in the southern lobe of the radio source. Section [4] 
deals with the gamma-ray production. Our results are 
there presented in the form of spectral energy distributions 
(SEDs), for different sets of parameters. We close with a 
brief discussion and a summary in Section [5j 

2. The protostar IRAS 16547-4247 and its host 
environment 

The source IRAS 16547—4247 corresponds to a young mas- 
sive star-forming region, associated with an O-type proto- 
star, located at 2.9 kpc (Garay et al. 2003). The luminosity 
of the source is L ~ 6.2 x 10 4 L© w 2.4 x 10 38 erg s _1 , 
peaking at the infrared, which makes of it the most lu- 
minous detected young stellar object (YSO) with thermal 
jets. Brooks et al. (2005) have reported the detection of a 
chain of H2 2.12 /im emission knots that trace the colli- 
mated outflow that emanates from the center of the source. 
The SED of the IR emission can be described by a modified 
blackbody function with a peak temperature of 30 K. The 
total mass of the cloud is M c \ = 9 x 10 2 M© (Garay et 
al. 2003). Molecular line observations indicate that the size 
of the cloud is ~ 0.38 pc in diameter (« 1.1 x 10 18 cm). 
If we assume a spherical geometry, the averaged particle 
(atoms of H) density of the cloud is n c \ « 5.2 x 10 5 cm~ 3 . 
The energy density of IR photons in the cloud, assuming an 
homogeneous distribution, is w p h w 1.8 x 10~ 9 erg cm~ 3 . 

The radio observations made by Garay et al. (2003) with 
the ATCA and the deeper observations by Rodriguez et al. 
(2005) with the VLA show the existence of a triple radio 
source inside the molecular cloud. The three components 
of the radio source are aligned in the northwest-southeast 
direction, with the outer lobes separated from the core by 
a projected distance of 0.14 pc. The central source is elon- 
gated and has a spectral index of 0.33 ± 0.05, consistent 
with free-free emission from a collimated jet (Rodriguez et 
al. 2005). The radio lobes have some substructure. The inte- 
grated emission from the northern lobe has a spectral index 
of — 0.32±0.29, of dubius thermal/non-thermal nature. The 
spectrum of the southern lobe radiation, instead, has an in- 
dex a = —0.59 ± 0.15. Actually, the shock heated material 




Fig. 1. Sketch of the scenario discussed in this work. 



can radiate enough as to ionize the surrounding medium 
and free-free absorption could modify to some extent the 
radio spectrum. However, we have no enough data to char- 
acterize the region in order to derive its thermal properties 
or the ionization degree. Otherwise, diffusive shock acceler- 
ation by strong non-relativistic shocks naturally produces 
a power-law particle distribution with an index similar to 
what is inferred from the observations. Thus, we take as a 
first order approximation the observed radio spectrum as 
the original one, which would correspond to a non-thermal 
particle distribution. The inferred linear size for this lobe 
is « 1.1 x 10 16 cm. Non-thermal radio emission has been 
associated in the past with the outflows of a few massive 
YSO (see Rodriguez et al. 2005, and references therein). In 
any case the molecular cloud is so massive and luminous as 
in the case of IRAS 16547-4247. The flux density of the 
southern lobe is 2.8 ± 0.1 mJy at 8.46 GHz. In Figure □ a 
sketch of the scenario discussed in this paper is shown. 

3. Particle acceleration and losses 

At the termination point of the thermal jet a strong shock 
front is expected to be formed. Charged particles can be 
accelerated at the shock, with an acceleration rate oc r]B 
(see below), where B is the magnetic field and the efficiency 
77 depends on the acceleration mechanism and its details. 
In the case of diffusive shock acceleration 77 ~ (w s /c) 2 // sc , 
where v s is the shock velocity and / sc is the ratio of the 
mean free path of particles to their gyro-radius. Close to 
the Bohm limit / sc ~ 1. The terminal velocity of the colli- 
mated outflow in IRAS 16547-4247 is unknown, but Marti, 
Rodriguez & Reipurth (1995) have determined a velocity in 
the range 600-1400 km s" 1 for the HH 80-81 thermal radio 
jet, which is powered also by a massive YSO. Adopting a 
value of - 1000 km s^ 1 for IRAS 16547-4247, we have a 
reasonable efficiency of r\ ~ 10~ 5 for the southern lobe of 
this source. The efficiency in the northern lobe would be 
lower, since the non-thermal emission does not seem to be 
dominant there. 

The average matter density of the cloud is 5.2 x 
10 5 cm -3 . It is noted nevertheless that in the shocked ma- 
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terial the density should be ~ 4 times this value. In our 
calculations of the Bremsstrahlung and pion-decay radia- 
tion we have adopted this higher density. 

As mentioned above, the rate of energy gain for elec- 
trons at the acceleration region is: 



7. 



r/eBc 



e, gain — q : 

m e c z 



(1) 



where 7 is the Lorentz factor of the particle, and m e is the 
electron rest mass. Similarly, for protons, 

- T]eBc o\ 

7p, gain — 9 • \£) 

m p c z 



3.1. Non-thermal radiation losses and maximum energies 

In order to get the maximum energy of the primary parti- 
cles we have to balance the energy gain and loss rates. In 
the case of electrons, the relevant losses are synchrotron, 
IC, and relativistic Bremsstrahlung losses. The latter two, 
under the conditions described in the previous section, are 
given by: 



Tic - -3.2 x Kr 17 7 2 s-\ 



7Br 



-7.4 x 1(T 10 7 s _1 . 



(3) 
(4) 



The expression for the synchrotron losses depends on the 
magnetic field B and is given below. We can estimate the 
magnetic field assuming equipartition in energy density be- 
tween field and particles in the radio lobe, so that: 



8vr 



(5) 



Here, u Cl is the energy density of primary electrons, u p that 
of primary protons, and it e2 corresponds to secondary pairs 
that will result from the charged pions produced in inelastic 
pp collisions. In each case (i = ei,p, 62): 



Ein(Ek)dEi, 



(6) 



where n(-E,) = KiE~ i exp(— _E i /£'- nax ) is the particle den- 
sity distribution [cm -3 erg ]. 

It is a fact that there are relativistic electrons or pairs 
in the southern lobe. This population of relativistic lep- 
tons can be produced via acceleration either of electrons 
or protons. In the former case, the accelerated electrons 
would be the responsible for the production of the detected 
synchrotron emission. In the latter case, accelerated pro- 
tons would suffer pp interactions with nuclei of the ambi- 
ent medium generating charged pions that would decay to 
muons, and these subsequently decaying to pairs. Then, 
these pairs would be the emitters of the observed non- 
thermal radiation. Although the calculation of the neutrino 
emission is out of the scope of this work, we note that neu- 
trino luminosities would be similar to those of 7r°-decay 
gamma-rays. 

The most likely scenario is one with acceleration of both 
electrons and protons. Nevertheless, the relative number of 
accelerated protons is unknown, so we will write u p — au Cl 
and then u C2 = fu p . The value of / was estimated using 
the average ratio of the number of secondary pairs to ir°- 
decay photons (Kelner, Aharonian & Bugayov, 2006). We 



Table 2. Break energy in the electron spectra for different 
cases discussed in the text. Normalization constants for the 
power-law distributions above the break energy are shown 
as well. 



a 


E h 
[eV] 


K' ei 
[erg 1 " 01 cm 


- 3 ] 


K' C2 
[erg r °2 C m 3 ] 





9.4 x 10 a 


3.7 x 1(T 


11 




1 


3.1 x 10 10 


9.8 x 1(T 


11 


3.2 x 10- 11 


100 


2.2 x 10 10 


9.5 x 10- 


13 


5.4 x 1Q- 11 



will consider three cases: a = (no proton acceleration), 
a = 1 (equal energy density in protons as in electrons), and 
a = 100 (proton dominance, as it is the case in the galactic 
cosmic rays). 

For a distribution of protons n(E p ), the energy spectra 
of secondary electron-positron pairs was calculated using 
the new parametrization of the inelastic cross-section of pp 
interactions given by Kelner et al. (2006): 



crinei(-Ep) = 34.3 + 1.88L + 0.25L 2 mb, 



(7) 



where L — \n(E p /l erg). In the cases a = and a = 1, 
r ci = 2.18, from the radio observations. For a = 1, when 
primary electrons are still dominant, the same slope is as- 
sumed for them and for protons. Then, for the secondaries, 
we get r c2 = 2.13. In the case of a = 100, where the sec- 
ondary leptonic emission dominates, r c2 must be 2.18 and 
then r p = r ci = 2.27 for the primary particles. 

Using the standard synchrotron formulae (Ginzburg 
& Syrovatskii, 1964) and the data from Rodriguez et al. 
(2005) , along with the equipartition condition and the equa- 
tions given above, we get the magnetic fields and the nor- 
malization constants shown in Table [TJ for each case con- 
sidered. Typically, B is about few lO -3 G. 

Now, with this magnetic field we can calculate the syn- 
chrotron losses in each case: 



1.0 x 1Q- 14 7 2 s~\ 



7synchr(a = 0) 
7synchr(a = 1) 
7synchr(a = 100) = -3.3 X 10 



1.1 x 1Q- 14 7 2 s~\ 



" 14 7 2 b- 1 . 



(8) 
(9) 
(10) 



The different losses for electrons are shown in Fig. [21 
along with the acceleration rate, for a value of B of 
2.5 x 10 -3 G (case a — 1). We can see that low energy 
electrons cool mainly through relativistic Bremsstrahlung, 
whereas the losses of the most energetic ones are dominated 
by synchrotron radiation. The dominant cooling regime 
changes at E e = E\> , where E e is the electron energy, either 
of primaries or secondaries. For a steady state, under the 
conditions considered here, the electrons will have an index 
r eii e 2 for energies E e < E\> and r oi!C2 + 1 for E e > E^. 
In Table [5] we show the values of E\, for the different cases 
under consideration, along with the corresponding normal- 
ization constants at high energies. 

The relevant losses for protons are from inelastic pp col- 
lisions with the cloud material. In our case, these losses are 
approximately given by (Mannheim & Schlickeiser 1994): 



Ipp 



-9.2 x 10" 



0.95 + 0.06 In 
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Table 1. Magnetic field and normalization constants for the different cases considered in the text. 
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Fig. 2. Energy loss rates for synchrotron and relativistic 
Bremsstrahlung processes. The IC energy loss rate is not 
shown since it is negligible. The acceleration rate and the 
particle rate of diffusion out of the emitting region are also 
shown. 



where now 7 is the Lorentz factor of the protons. Equating 
with Eq. ((2]), we get the maximum energy of the accelerated 
protons, which results of the order of ~ 10 15 eV. Although 
the gyro-radius of these particles can be contained in the 
acceleration region, they can diffuse out of the cloud be- 
fore reaching the mentioned maximum energy, hence a more 
tight constraint on the maximum energy can be obtained 
from the following equation: 



Table 3. Maximum energies obtained for electrons and pro- 
tons accelerated in the southern lobe of the radio source 
associated with IRAS 16547—4247. The maximum energy 
for secondary pairs is shown in the last column. 
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context, the spectrum of the particles is determined by the 
following transport equation (e.g. Khangulyan et al. 2007): 

dn {t, 7) jdt + djn (t, 7) /d 1 + n (t, 7) /r esc = Q(t, 7), (13) 

where Q(t, 7) is the function for the particle injection, which 
takes place during the lifetime of the source (ruf e ), esti- 
mated in > 100 yr« 3 x 10 9 s (see Garay et al. 2003). Here, t 
is the time, and 7 is the lepton Lorentz factor. Admittedly, 
the age of the source is not well constrained, although it 
could hardly be significantly younger than ~ 100 yr, the 
minimum jet crossing time, since the shock has been prob- 
ably active for most of this period while the shocked mate- 
rial was being displaced to its present location. Concerning 
r osc , at this stage we take it as being the shortest diffusion 
timescale, i.e. that corresponding to the maximum energy 
particles. The time-derivative 7 is a function accounting 
for all the energy losses affecting leptons, i.e. basically syn- 
chrotron and relativistic Bremsstrahlung. The solution of 
Eq. HID is: 



r/eBc 



lobe 

D 



(12) 



where D is the diffusion coefficient and ri o b w 1.1 x 10 16 
cm the size of the acceleration region. This assumption is 
supported by the fact that particles lose all their energy 
before they diffuse out of the emitting region (i.e. T coo i ing < 
Tdiff) (see Fig. [2]). Working on the Bohm limit we can set 
D = Db = f g c/3, with r g the gyro-radius. In such a case we 
get the maximum energies shown in Tabled which are more 
than one order of magnitude lower than in the previous 
calculation. 



3.2. Relativistic particle distribution 

The lepton particle distribution in the emitting region can 
be calculated, as a first order approximation, adopting a 
one-zone model, in which particles, once injected following a 
power-law energy distribution, evolve suffering mainly syn- 
chrotron and relativistic Bremsstrahlung losses. In such a 



i(«,7) = \ J Q{t-r,i)e- r My^di, 
7 



(14) 



where 

7eff 

t = 



and r(7,7') = 
7' J 7 

7 7 



dj' 
11 



(15) 



The time dependence comes through Eq. (fT5|) . For small t, 
7 7cff ; for larger t, 7 <C j c s- Since y e g can be arbitrarily 
large, for 7' in Eq. (| 14[) above the maximum Lorentz factor 
Q will be zero. In addition, if r, being shorter than r coo ii ng , 
becomes larger than r osc , the final particle spectrum will be 
affected by particle escape. 

The computed particle energy distributions of primary 
and secondary leptons are shown in Fig. [3J As seen in the 
figure, for injection timescales T ln j ~ 10 9 s< ri;f e , the par- 
ticle distribution n (t, 7) has already reached the steady 
regime, as expected from the fact that r coo ii ng ~ 10 9 s. 
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Fig. 3. Spectral energy distributions of primary electrons and secondary pairs. Distributions for different particle injection 
timescales (ri n j) are shown in each panel. Different curves correspond to the following values of r; n j: 10 7 s (long-dashed 
line, green), 10 s s (dot-dashed line, blue), 10 9 s (dotted line, red) and 10 10 s (solid line, black). Notice that for ages 
> 10 9 s the steady regime is reached (the curves appear in color only in the on-line version of the paper). 



Moreover, the final particle energy distribution is not af- 
fected by r csc (which is actually a lower limit for the 
timescale of diffusion particle escape), since r csc Jj r coo ii ng 
(see Fig. [2]), i.e. particles will radiate inside the emitting 
region. 

Regarding protons, their diffusion timescales and source 
lifetime are long enough to let them radiate almost all their 
energy in the hot spot. Moreover, the energy loss timescale 
is only slightly dependent on energy. This implies that pro- 
ton energy distribution, in the context of our scenario, keeps 
almost the same spectral shape as that of the injected one, 
and reaches, as it is the case for leptons, the steady regime, 
without suffering significant impact from escape losses. 

3.3. Thermal radiation and ionization losses in the shocked 
region 

In the context of an outflow interacting with its environ- 
ment, a strong shock of several hundreds km s _1 can heat 
the shocked cloud material up to temperatures of ~ 10 7 K, 
producing emission that would peak at soft X-rays. The 
heating of the medium would take place via shock compres- 
sion, and ionization and Coulomb scattering of the thermal 
atoms/ions would occur due to the presence of relativistic 
leptons and protons. These soft X-rays should be strongly 
absorbed by the surrounding medium, provided the large 
hydrogen column density of the cloud. Hence, this radiation 
will be hard to detect. In any case, this would lead to ion- 
ization of the cold material surrounding and/or embedded 



by the emitting region, which may affect to some extent 
the radio emission via free-free absorption. Since the ob- 
served radio spectrum seems to be non-thermal and we aim 
at obtaining just hrst order estimates of the high-energy 
emission, we assume here that the overall effect of the free- 
free absorption is not dominant in the southern lobe of 
the source. Nevertheless, a more detailed characterization 
of the original radio spectrum would deserve further study. 
Another observational consequence of the ionization of the 
medium would be the detection of recombination lines, al- 
though their study is beyond of the scope of our work. 

It is worth mentioning that, if the emitting region were 
not significantly ionized by the jet termination shock, ion- 
ization losses would affect the particle energy distribution. 
The effect of this loss channel, which energy loss rate does 
not depend on energy, would be to harden the electron en- 
ergy distribution at the energies where ionization losses are 
greater than those due to relativistic Bremsstrahlung, i.e. 
at Lorentz factors 7 < 10 3 . We note that, if energy losses 
were indeed dominant at 7 < 10 3 , their impact would be 
a pretty hard radio spectrum, which is apparently not the 
case. 



4. Production of gamma-rays and lower energy 
radiation 

From the previous section it is clear that at the termination 
shock there will be particles energetic enough as to produce 
gamma-ray emission. In this section we will use the particle 
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Fig. 4. Spectral energy distribution for a pure leptonic case. 

distributions, ambient photon density and magnetic fields 
already estimated to calculate the non-thermal spectral en- 
ergy distributions expected from the different particle inter- 
actions in the southern lobe of the radio source associated 
with IRAS 16547-4247. 

4.1. Leptonic interactions 

Synchrotron, rclativistic Bremsstrahlung and IC emissivi- 
ties for primary and secondary leptons were calculated us- 
ing the standard formulae (e.g. Blumenthal & Gould 1970, 
Pacholczyk 1970). Full Klein- Nishina cross section was used 
for the IC computation. The emission region was considered 
spherical, with a radius r = Hobo- FigureHJshows the results 
for the pure leptonic case (a = 0). 

Focusing on the energy bands at which non-thermal 
emission may be detectable, we see that the high-energy 
emission is dominated by the relativistic Bremsstrahlung, 
with a peak of ~ 10 32 erg s _1 at E 1 ~ 1 MeV. At energies 
E 1 Jb 1 GeV, the source presents luminosities of ~ 10 31 
erg s _1 , with a cutoff ~ 10 GeV. At X-rays, in the range 
1-10 keV, the expected luminosities are about 10 30 erg s _1 , 
presenting a softening of the spectrum due to the high en- 
ergy cutoff of the electron energy distribution. 



4.2. Hadronic interactions 

In the cases with a = 1 and a = 100, protons will have 
a significant contribution. The differential 7-ray emissivity 
generated through 7r°-decay was calculated as 



q 7 {E 7 ) = 2 



<2V° 



where E™ in (E 7 ) = E 1 



cIEt, 



(16) 



Applying the 5-function 



j 1) — -r 4E _^ 

approximation for the differential cross section^ Aharonian 
& Atoyan, 2000), the pion emissivity becomes 

q 7r o{E- jr ) = 47r Jj^ 5{E„ - KE Mn )J p (E p ) a pp (E p ) dE p 



1 This approximation considers only the most energetic neu- 
tral pion that is produced in the pp reaction aside of a fireball 
composed by a certain number of less energetic 7r-mesons of each 
flavor. See a discussion in Pfrommer and En/31in (2004). 
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Fig. 5. Spectral energy distribution for a mixed primary 
population of relativistic electrons and protons. We indicate 
with numbers (1) and (2) the contribution from primary 
and secondary electrons, respectively. 

= f J P K c2 + ^) ^ K^ 2 + sc) (17) 

for proton energies greater than the energy threshold E t h — 
1.22 GeV and lower than 100 GeV. Here, k is the mean 
fraction of the kinetic energy -Ekin = E p — m p c 2 of the 
proton transferred to a secondary meson per collision. For 
a broad energy region (GeV to TeV) k ~ 0.17. The total 
cross section of the inelastic pp collisions is given by Eq. ([7]) . 
In Eq. mi) J p {E p ) = (Air / c)n(E p ) is the proton flux. For 
E p > 100 GeV, the equations given by Kelner et al (2006) 
were used. 

The specific luminosity can then be estimated as: 



EjLe^ 



E~l \ g 7 (E 7 )n(r)dr, 



(18) 



where n(r) is the ambient particle density. 

Results from our calculations can be seen in Figures [5] 
andO The hadronic emission is similar in luminosity to the 
relativistic Bremsstrahlung (~ 10 32 erg s _1 ), but extends 
up to higher energies ~ 1 TeV. Primary electrons dominate 
the synchrotron and relativistic Bremsstrahlung emission 
in the case a — 1, but for a = 100 (see Fig. 4) the sec- 
ondary leptons produce the most important contribution. 
In neither case the IC plays a significant role. Basically, 
the emission above ~ 1 GeV is due to neutral pion de- 
cays, whereas the X-ray and soft gamma-rays are produced 
by relativistic Bremsstrahlung. For the cases a — 1 and 
100, the synchrotron and relativistic Bremsstrahlung com- 
ponents produced by secondary pairs are similar in spectral 
shape to those generated by primary electrons in the same 
cases, but with different behavior at low energies, due to 
the low-energy cutoff introduced by the minimum energy 
at which secondaries are injected. 



5. Discussion and summary 

Gamma-ray sources on the galactic plane are usually as- 
sociated with star forming regions (Romero, Benaglia & 
Torres 1999). However, the sources detected by EGRET are 
far more luminous than the sources that can be produced 
by massive YSOs according to our calculations. EGRET 
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Fig. 6. Idem Fig.O but for a case dominated by relativistic 
protons. 



the region, as well as the relatively high kinetic luminos- 
ity of the jet, render likely that radiation will be signif- 
icantly generated by channels other than synchrotron ra- 
diation, like relativistic Bremsstrahlung or even pp inter- 
actions. After studying in a first-order approach a reason- 
able scenario for the non-thermal emitting region, we con- 
clude that YSOs like IRAS 16547-4247 could produce de- 
tectable non-thermal emission via synchrotron, relativistic 
Bremsstrahlung and pp interactions from radio up to very 
high-energies. This opens a new window for studies of mas- 
sive star formation in the near future. 
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unidentified sources are most likely produced by pulsars, 
supernova remnants, massive stars (Romero 2001), and mi- 
croquasars (Bosch- Ramon, Romero & Paredes 2005). The 
emission of massive YSOs falls below the sensitivity of 
EGRET, but they can be detected in principle by GLAST. 
At very high energies (E^ <; 100 GeV), Cherenkov ar- 
rays like HESS-II could be able to detect a source like 
IRAS 16547—4247. Since such instruments might find the 
high-energy cutoff, valuable information could be obtained 
about the true efficiency of the acceleration process. 

At X-rays, the source may be detectable by Chandra 
and XMM-Newton, probably as point-like, through a deep 
observation of the region. Briefly looking at archive XMM- 
Newton 30 ks observations of IRAS 16547-4247, there is 
not significant emission above the background from the 
southern lobe region. This yields an upper-limit for the 1— 
10 keV luminosity of ~ 10 31 erg s _1 at the distance of 
2.9 kpc. This value is not in contradiction with the lumi- 
nosities of synchrotron and the relativistic Bremsstrahlung 
components expected in this energy range (see Figs. HI [5] 
and [6]) . Nevertheless, a deeper observation may detect ra- 
diation from this region. In case the temperatures of the 
shocked material were high enough, a thermal component 
may be detectable in the X-rays as well. In the UV-band, 
the luminosity may still be higher than that produced by 
non-thermal particles, but could be difficult to detect. 

An issue that is relevant for the feasibility of the sce- 
nario discussed in this paper is that the non-thermal lu- 
minosities are well below the kinetic energy budget of the 
jet in IRAS 16547-4247. The jet mass-loss rate has been 
estimated by Garay et al. (2003) in ~ 10~ 5 M yr" 1 . With 
a velocity of ~ 10 3 km s _1 , the jet kinetic luminosity is 
~ 3 x 10 36 erg s^ 1 , more than three orders of magnitude 
larger than the non-thermal luminosity. Moreover, the en- 
ergy density of the shocked thermal gas should be about 
two orders of magnitude larger than that of the non-thermal 
particles. Therefore, to sustain the non-thermal emission es- 
timated in this work under reasonable assumptions seems 
not to be a problem from the energetic point of view. 

Summarizing, the fact that synchrotron radio emission 
has been detected from IRAS 16547—4247 already indi- 
cates that there is conversion of energy of some sort to 
non-thermal particles, and the high densities present in 
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